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1. SUMMARY 


Mathematical modeling for the description of the origin 
and disposition of combustion-generated pollutants in gas 
turbines is presented. A unified model in modular form is 
proposed which includes kinetics, recirculation, turbulent 
mixing, multiphase flow effects, swirl and secondary air 
injection. Subelements of the overall model have been applied 
to data relevant to laboratory reactors and practical combustor 
configurations. Comparisons between the theory and available 
data show excellent agreement for basic CO/F^/Air chemical 
systems. For hydrocarbons the trends are predicted well 
including higher- than-equilibrium NO levels within the fuel - 
rich regime. Although the need for improved accuracy in fuel- 
rich combustion is indicated, comparisons with actual jet- 
engine data in terms of the effect of combustor-inlet temperature 
is excellent. In addition, excellent agreement with data is 
obtained regarding reduced NO emissions with water droplet and 
steam injection. 

Comparisons of predictions made on swirling flame data 
show that while good agreement can be obtained additional work 
is required in terms of the turbulent transport properties. 

In general, the modeling presented herein has aided in the 
understanding of the emissions problem and has delineated areas 
needing further study. 



2. INTRODUCTION 


Recent studies ^ ^ ^ ) have established that the density of 
pollutant emissions from aircraft turbine engines, particularly 
in the vicinity of airports, is a matter of increasing public 
concern and will become even more important in future years. 

This is a problem which has stimulated the output of a large 
body of data on the emission characteristics of existing as 
well as advanced combustor concepts (4,5, 6, 7,8) ^ However, the 
data itself has not been sufficient to explain all the trends 
exhibited by the observations. Part of this difficulty is due 
to instrumentation problems and part has been due to the lack 
of appropriate mathematical models useful in providing the 
insight necessary to understand the data. The emissions problem 
arises out of past engine developmental procedures which have 
emphasized overall combustion efficiency and performance while 
giving little attention to the resulting pollutant generation 
and emissions. 

The conventional treatment of the combustor process 
involves one-dimensional concepts and only recently has some 
consideration been given to the inclusion of the kinetics 
process as an improvement over the equilibrium assumption. 

Very little work has been done in the past to couple valid 
kinetics with the appropriate fluid mechanics. 

It is generally recognized that the pollutant emission 
characteristics of aircraft turbine engines are related to the 
failure of the exhaust gases to attain chemical equilibrium. 
Analysis of exhaust gas samples indicate that the departure from 
equilibrium, especially for species such as nitric oxide, CO and 
hydrocarbons, is extremely high. In actual practice, the effect 
is to produce greater quantities of these pollutants than would 
exist if equilibrium prevailed throughout the engine. Thus, the 
exhaust composition of these pollutants is dependent upon the 
rate limiting processes within the combustor. In order to formu- 
late a realistic mathematical model of these processes, a knowl- 
edge of the detailed flow field within the combustor is required. 

The combustion process in such annular and canannular combus- 
tors can be described in terms of the flow pattern in two zones; 
a primary and secondary, which are shown schematically in Figure 1 
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In the primary zone, air enters through swirlers (bladed 
passages which impart a swirl to the air flowing through) 
around the fuel spray nozzles and air holes in the combustor 
liner. In this zone, the fuel and air are mixed and partially 
burned at near stoichiometric conditions ( $> = 0.8 - 1.2 ). 
However, the discrete injection of the fuel results in a wide 
spectrum of equivalence ratios in this region. This provides 
for nonuniform fuel evaporation and mixing giving rise to 
locally fuel-rich and fuel-lean domains. This equivalence 
ratio variation together with the associated temperature 
variation will lead to: (1) partial oxidation and thermal 

cracking (pyrolysis) with the attendant formation of hydro- 
carbon fragments, CO, and soot in the fuel-rich regions, 

(2) maximum formation of N0 X on the lean side, near stoichio- 
metric regions, and (3) quenching in the colder very lean 
regions. In practice, the concentra tion levels of these species 
depends upon the mixing and kinetic rates relative to the local 
residence times. A mathematical model representing the overall 
primary zone must account for the mechanisms controlling these 
rates as a function of initial fuel and air properties. 

In the secondary zone. Figure 1, the main bulk of the air 
is added and the reaction is continued. Further dilution occurs 
in the secondary zone to reduce the product gas temperature to 
the desired turbine inlet conditions. For an overall combustor 
equivalence ratio of about 0.25, approximately 75 percent of the 
total air flow is introduced in the secondary zone. Part of the 
air from the secondary holes recirculates upstream into the 
primary zone and improves the combustion stability and part is 
injected close to the flame-tube wall providing film cooling 
(Figure 1). 

The assumption of one-dimensional flow is inadequate in the 
treatment of the secondary zone because air is being injected 
through holes and cooling slots along the combustor length. 

The one-dimensional flow assumption which implies flow properties 
are uniform across any cross-section requires the very poor 
hypothesis that the injected air mixes instantaneously with the 
bulk flow. In reality, the injected air can require appreciable 
distances to mix with the bulk flow. Thus, in addition to the 
nonuniform initial condition generated by the primary zone flow 
field, the boundary conditions along the secondary zone walls 
will force a persistence of the nonuniformity throughout the 
chamber. The effect of this nonuniformity includes possible 
quenching of the reactions due to the cooling near the walls. 
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Again, as in the primary zone, the velocity, concentration, and 
thermal fields in the secondary zone depend upon the relative 
importance of the mixing, chemical and flow times. Thus, a 
mathematical model for the secondary zone must include the 
inter-relationship of these processes and take account of 
their dependence upon the primary zone characteristics as 
well as the dilution and coolant air distributions and states 
of flow. 

Thus, the objective of this work was to develop a modular 
computer program for the description of gas turbine combustor 
flow and chemical fields for the prediction of combustion 
efficiency and pollutant formation and disposition. 


3. MATHEMATICAL MODEL FORMULATION 


a. General Model Requirements 

The jet engine combustor involves perhaps one of the most 
complex reacting flow fields of practical interest. The com- 
plexity arises out of the number of discrete mechanisms which 
occur simultaneously including heterogeneous reaction kinetics 
and turbulent mixing with recirculation patterns all taking 
place in a subsonic flow. Such a problem is in general elliptic 
and would require a Navier-S tokes treatment to thoroughly repre- 
sent the flow field. To completely model the complex flow 
pattern existing in the primary zone of a real combustor is 
beyond the capability of present analytical methods. However, 
several approximate models can be postulated which character- 
ize this zone in terms of mixing rates, reaction rates, and 
residence times. For example, one limit of operation of this 
zone is characterized by a "stirred reactor." That is, a well 
established recirculation zone exists and the process in terms 
of the degree of reaction is controlled by the chemical kinetic 
rates. Another limiting representation of this zone is char- 
acterized by a discrete jet mixing and combustion process. In 
this case, a parabolic analysis would be appropriate. The more 
general representation involves a combination of stirred reactor 
and jet type flows as shown in Figure 1. That is, a local region 
involving a recirculation, or stirred reactor zone, surrounded 
by a jet type mixing and combustion zone. Fuel injection and 
spray patterns must be accounted for in terms of the nozzle 
characteristics. 
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The analysis of the secondary zone must include the 
following principal effects: " 

1. Turbulent mixing and diffusion of the gases, 
droplets (if they should still exist), and 
carbon particles, including the effect of 
mixing of the injected air (injected through 
the cooling slots and dilution holes) and the 
bulk reacting gases; 

2. Finite-rate chemical kinetics which describe 
not only the reactions between the combustion 
products (CO 2 , CO, N0 X , etc.) but also the 
further oxidation of the unburned hydrocarbon 
and of the carbon particles. 

b. The Modular Concept 

The modular concept proposed here includes a rather general 
set of discrete mechanisms which control combustion efficiency 
and pollutant emissions. These mechanisms can be conveniently 
represented in terms of chemical kinetics and fluid mechanics 
processes as shown in Figure 2 . However, it should be noted 
that these ingredients are coupled in the actual combustor. 

Thus Figure 2 shows that the reaction kinetics include droplet 
combustion and involves the exothermic reactions associated with 
the production of CO 2 and H 2 O, the non-exo thermic reactions 
involving combustion intermediates, as well as N0 X and S0 X 
(appropriate for sulfur bearing fuels), and finally, CO, C x Hy , 
and soot, C( s ) , which can be associated with the combustion 
efficiency. In terms of pollutants, the net output of the 
chemical reaction processes are the concentrations of CO , 

C x Hy , C( s ), N0 X and S0 X . Of course, the reaction processes 
are ocurring with a complex flow field involving fuel droplets 
(and soot) suspended in the gas. Thus the multiphase fluid 
mechanics includes both gas and droplet mixing in the presence 
of swirl and recirculation. 

The primary virtue of the proposed model lies in its 
inherent ability to respond to changes in many of the operating 
conditions and geometric factors which are of practical interest. 
This capability is summarized in TABLE I. 



In general, the model will predict the details of the flow 
including the local velocity components, concentrations of gas 
and droplets, temperature and pressure (for a specified liner 
contour). In addition, mean properties within a recirculation 
zone are predicted. These variables are displayed in Column I 
while Column II shows their initial conditions which are 
arbitrarily specified. Columns III, IV and V detail the boundary 
conditions where it is worth noting that the area contour or 
wall pressure distributions may be specified. This option is 
of particular interest for advanced high speed combustors where 
sufficient dynamic head is available for significant pressure 
conversion. Thus, it may be desirable to determine a liner 
contour for some specified pressure distribution that provides 
control of the velocity field as well as of the kinetics. In 
particular, control of mixing rates, residence times, and com- 
bustion rates can be studied with the option as outlined in 
Column III. Column IV outlines the treatment of slots and holes. 
The "lateral" injection process through the holes is treated as 
a unit problem for the determination of the penetration into 
the main flow at the particular combustor cross section. The 
depth of penetration defines the radial location at which the 
injected flow becomes essentially parallel to the main com- 
bustor flow. Here, the flow is circumferentially mass averaged 
after taking into account the acceleration of the main combus- 
tor flow due to the area occupied by the injected air. The 
"new" profiles of the flow variables provide the information 
required to continue the solution downstream. 


c. Chemical Kinetics 

It is well established that in high performance systems 
the chemical equilibrium assumption is not valid for most species 
of interest. This is particularly true of N0 X even for residence 
times well beyond the millisecond range. Carbon monoxide is a 
problem since in many systems its oxidation to CC >2 is halted due 
to quenching at reduced temperatures or dissociation at peak 
temperatures. Only a few points of combustion inefficiency are 
associated with unacceptable CO emission levels. This is also 
true of the unburned hydrocarbons (UHC). 

One of the most significant findings associated with 
recent pollution-oriented kinetics modeling is the extreme 
importance of many of the intermediate, atoms and free radicals 
that are an inherent part of any oxidation process. Reference 9. Super- 
ecruilibrium levels of such soecies occur during the early stages 
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of combustion and they persist, due to the relatively slow three - 
body recombination reactions/even beyond what is normally con- 
sidered the overall reaction time. This behavior accounts, at 
least in part, for such observations as "Prompt NO", Reference 10. 

Perhaps the single most important consideration in the formu- 
lation of appropriate kinetic mechanisms is the recognition of the 
spectrum of local states that are encountered in various combustor 
flow fields. This includes the distribution of fuel/air ratios, 
temperatures, pressure ratios, and residence times. A rather 
general consideration must be given to species and reactions 
each of which may become important under the varying conditions 
encountered in practical combustor configurations. 

It is common to distinguish between combustion reactions 
involving the exothermic formation of CO 2 and H 2 O , and the NOx 
formation reactions. In practice, however, the two processes 
cannot be distinguished and they are coupled. For example, 
under many combustion conditions of interest the Zeldovich 
mechanism appears to govern the formation of NO through the 
rate controlling reactions: 


O + N 2 ^ NO + N 


( 1 ) 


Now, the assumption that O is at its equilibrium level is 
generally not valid. Superequilibrium levels of atomic oxygen 
as well as other intermediates are predicted which is totally 
consistent with observed ignition delay and reaction time data. 
In view of these observations we have developed models which 
treat the combustion and N0 X formation mechanisms within the 
framework of a unified scheme. All reactions are coupled and 
neither equilibrium nor steady- state assumptions are made. 

Our early work. Reference 11, demonstrated that many 
observations on the rate of combustion of hydrocarbons could be 
predicted by a relatively simple kinetic scheme. We termed the 
mechanism "quasi-global " which has as a key element a subglobal 
partial oxidation step: 


C H 
n m 




+ nCO 


( 2 ) 


This reaction is unidirectional with an empirically determined 
rate (grams of fuel /cc/sec) given by 

5.52 (10) 8 (p) - * 825 (T) [C H l* 5 [0_] exp (-12,200/T) (3) 

. n m Z 

where p must be given in- atmospheres , T in degrees Kelvin and 
[ ] denotes molar concentration. 
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Coupled to this subglobal step are the intermediate reversible 
reactions given in TABLE II. Some observations on this scheme 
should be noted. First, it is generally agreed that in hydrogen 
bearing systems carbon monoxide is most rapidly oxidized to CO 2 
via the reaction 

CO + OH ~^..C0 2 + H (4) 

We have, however, included other reactions involving CO and 
oxygen for completeness including, 

(5) 

( 6 ) 

where M is the general third body. Reactions 5 and 6 are much 
slower than reaction 4 but their inclusion was necessary for 
basic studies performed on systems including the CO/Air system 
such as reported upon in Reference 12. In addition, a number 
of reactions involving N0 X are included which represent a 
necessary extension of the basic Zeldovich mechanism to account 
for certain of the ambient long time N0-to-N0 2 conversion 
reactions which occur in the 'atmosphere particularly when 
coupled with appropriate daylight photochemical mechanisms. 

There are, however, other species of interest which can 
play an important role in the hydrocarbon oxidation process. 

A problem remains and this is mostly in fuel-rich hydrocarbon 
combustion. Hydrocarbon fragments, radicals, and a variety of 
oxygenated species are produced during the combustion process 
and accounting for them is of potential importance in adquately 
representing the fuel oxidation process as well as the N 0 X 
formation process. 

Perhaps one of the most exhaustively studied hydrocarbon 
oxidation mechanisms is associated with the methane/air system. 
References 13, 14 and 15. Although there is an apparent 

controversy over the importance of certain species and reactions 
governing the fuel oxidation process agreement on several aspects 
are worth noting: (1) A basic H 2 /Air mechanism is crucial, 

(2) CO oxidation is controlled by the hydroxyl radical (OH) con- 
centration level, and (3) the methyl radical (CH 3 ) dominates 
the initial phases of the oxidation process. The points in 
question relate to ( 1 ) the relevance of formaldyhyde , ( 2 ) the 

importance of hydrogen peroxide (H 2 O 2 ) and the hydroperoxyl 


CO + 0, 


co 2 + 0 


and 


CO + 0 + M^==^C0 2 + M 
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radical (HC> 2 ) ., and (3 ) the appearance of higher hydrocarbons 
such as ethane (C 2 H 5 ) as experimentally observed in Reference 13. 
Finally, little work: has been done to quantify the effects of 
other intermediates including CN and HCN type species in terms 
of. their influence on NOjj formation. Fenimore, Reference 10, 
has suggested that such species may play a role in the "early" 
formation of NO. 

To answer some of these questions we have initiated a 
study on certain aspects of the methane oxidation process. 

To date, the investigation has been limited to a scheme we 
have postulated and the scheme proposed by Bowman and Seery, 
Reference 14. These are given in TABLES III and IV. 

respectively. We have kept a number of reactions not included 
in Reference 14 but the major difference is our retention of 
formaldehyde : (HCHO) . As cited above, a question regarding the 
importance of the explicit appearance of HCHO exists and has 
led to deleting it by employing the subglobal steps. Reference 14. 



CH 3 + 0 2 = H 2 0 + CHO , 

(7) 

and 


CH 3 +0 = CHO + H 2 , 

(8) 

in place 

of the sets, - 



' CH 3 + 0 = HCHO + H , 

(9a) 


HCHO + H = CHO + H 2 , ’ 

(9b) 

and 


CH 3 + 0 2 = HCHO + OH , 

( 10 a) 


HCHO + OH = CHO + -.HO 

( 10 b) 

By 

retaining reactions 9 and 10, 

together with addi 


reactions involving HCHO, TABLE III not only provides a means 
of assessing its importance as an intermediate but also permits 
the prediction of the HCHO emission level which is an important 
pollutant entering into the atmospheric chemistry problem. 

Now, the above mechanisms have been investigated in a 
variety of flows comprising the elements of our modular concept. 
Emphasis here is on the "qua si-global " model. 
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Recirculation Zone Model 


The existence of recirculating regions in the. primary zone 
of turbine combustors provides a mechanism for rapid mixing. 

If the backrriixing is sufficiently intense this local zone will 
tend toward a homogeneous state. As the limit of essentially 
complete mixing is approached the principal mechanism defining 
the state of flow becomes the chemical kinetics of the fuel 
oxidation process. The residence time becomes the controlling 
parameter for a given fuel and air state entering such a recircu- 
lating zone . 


The equations describing the process are straightforward 
and are deduced directly by application of the conservation of 
mass, energy and species. The conditions are: steady state, 
homogeneity within the reacting volume, and that the efflux state 
is identical to the state existing within the reaction zone. 

The resulting describing equations are: 

Global Conservation of Mass: 


m = 


constant 




( 11 ) 


where k is the kth component and includes i gas-phase species, 
and j droplet types. A droplet type (or class) is defined by 
its composition (fuel, water, etc.) and its size. The super- 
scripts I and o refer to inflow and outflow, respectively. 

Conservation of Energy : 

h = h 1 + Q/m (12) 

where Q is the net rate of external heat addition to the 
reactor and the inflow of enthalpy is written to allow each 
component to enter the reactor with an arbitrary temperature. 
In particular, , 

h = £ a.h. (T) + £ a. h . (T. ) , (13) 

i 1 1 j 3 3 3 

£ cgh. {T 1 )+ £ a 1 h . (T 1 ) (14) 

, i i i . • J 1 3 


and 


li = 
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where 


T = Local gas-phase temperature 
th 

T_..= Local j droplet temperature 

til 

Injection temperature of the i gas-phase specie 

til 

Injection temperature of the j droplet 

To complete the specification of the thermal field some con- 
straint governing the conservation of droplet energy is required 
For the present it will be assumed that each droplet type is at 
its respective boiling-point temperature, i.e., T. = T. 



Conservation of Species: 


The conservation of species k requires that its net out- 
flow be equal to the rate of production due to chemical and 
phase transition processes. Formally, this principle gives 
k equations of the form: 


\ = °i 


V r *G 
+ T [W + 
m k 




( 15 ) 


where superscripts G and P refer to homogeneous gas-phase and 
dropwise production rates, respectively. 


Homogeneous Gas Phase Reactions 
The volumetric production rates are given as follows: 


R 


W. = W. £ ( v - v. )k^ p 


m N | or. i ip 
P H P| — 1 


P=1 


ip~Vf^ W. 

^ p 1 = 1 * x 


N 


1 - 


P N / a. J^ip ^if) 


where 


N N 

m = £ v ; ■ N = £ (v". -v! ) 


P ± ip 


P i iP ip 


c,p 1=1 ‘ 1 


( 16 ) 


( 17 ) 
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for a chemical system containing N gas-phase species entering 
into R elementary reversible reactions given by: 


N : 

L v ! M . _ 

i=i ip 1 nr 


N 

T V M . , p = 1 , ... R 

i=l 1P 1 


(18) 


Dropwise Combustion and Evaporation 


For now, the droplet types will be specialized to fuel 
and water. Two limiting modes of consumption will be considered 
for the fuel droplet: (1) diffusion-controlled combustion and 
(2) evaporation only. In both cases a detailed homogeneous gas 

phase kinetics mechanism will be in effect. Of course, for the 
water only evaporation is relevant. 

Fuel : 

The two reactions considered are: 

1. C H + (n + j)' 0_ nCO„ + “ H 0 (19) 

n m 4 2 2 2 2 


2 . (C H ) _ . . . 

n m Liquid 


(C H ) 
n m gas 


( 20 ) 


In both cases, the rate of consumption of the liquid fuel is 
given by diffusion-controlled theory; viz: 


<”c H )P - 
n m 


- 3/2 p d, 


C H 
n m 


a 


I ■ 


C H i 
n m 


C H 
n m 


( 21 ) 


where 


a. 


C H 
n m 


= d' 


C H l 
n m 


C H 
n m 


a 


C H l 
n m 


( 22 ) 


and d is the droplet diameter. In addition. 


X c H ■ - 8 ( C-6>C H (1+B C H » 
n m p n m n m 


(23) 
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where 


A (Y ) ■- 

O. g N C H 
2 3 'nr 


_ _B.P. \ 
T-T ! 
CHI 
n m ■ 


C H 
n m 


C H 
n m 


and k and Cp are the thermal conductivity and specific heat of 
the vapor, respectively; H is the heat of combustion, N is the 
stoichiometric 0/F ratio, 6 is the droplet density, Q is the 
latent heat of vaporization, and A is a factor which differ- 
entiates between diffusion-controlled combustion (A=l ) and droplet 
evaporation (A=0) . 

Now, for A = 1 we have: 

P (n + 4 ] W °2 P 

V ~ nW + mW ^ W C H ^ 

2 C H n m 


(W ) = 

°2 


2 * P 

£ (w ) F 

nW + mW rT C H 
C H n m 


<W H 0) E 
2 g 


IE TT 

2 -- H 2° (w ) P - 
nW o + mW v C H 
C H n m 


<v> p < 27) 


and for A = 0: 


(w c H ) P 
n m g 


- (W ) 
' C H 
n m 


In addition. 


(Y ) = z'cT , i = all gases (2 9) 

2 9 i i 

where (Yo 2 )g is the mass fraction of O 2 in the gas-phase subsystem. 
Water : 

Here the reaction is: 

(HO) — ► (HO) (30) 

2 “i ^9 


where P U 1 \ 

<v ;- 3/2p v ? 



and 


Vo " 8 f«> n (1+B H 0> 

2 p H 2 0 2 


(32) 


with 


B. 


h 2 ° 


V <t - t v > 

Q 

H 2° 


(33) 


and 


°H. 


d 3 

h 2 ° 


= d 


I 3 n 2° l 

H o 0 I 

2 V* 


(34) 


S ta te : 


P = 


rt L 

i W i 


, i = all gases 


( 35 ) 


where p is the mixture density and T is the gas-phase temperature 

Although there exist more or less standard techniques for 
the solution of such algebraic equations, an alternative 
approach is that of seeking the asymptotic (steady state) 
solution to a transient problem wherein the "boundary conditions' 
are held fixed at the desired steady state values. In essence, 
our approach involves the use of the non- steady form of the 
species conservation equations, viz: 


ffk 

dt 




(36) 


where t is the time variable of significance only during the 
transient period. Thus, Equation (36) is identical to 
Equation (15) when: 



do^ 

dt 


- 0 . 


(37) 
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The remaining working equations retain their steady state 
form. TO initiate a calculation requires the specification 
of the o^'s at t = 0. We start with the equilibrium state 
although this "initial" composition may be chosen arbitrarily. 

The above formulation represents only part of the overall 
combustor model. In general the recirculation zone is imbedded 
within a "directed" flow in the primary zone and the two zones 
are coupled through diffusive transport across the dividing 
streamline. Furthermore, in the downstream secondary zone the 
mixing and kinetics processes are continued and in some cases 
must include the effects of secondary air injection for cooling 
and dilution. 


e. Mixing Model 

The principal element forming the framework of the 
unified combustor model is the mixing analysis. It is within 
this framework that recirculation zones, additional kinetics, 
secondary injection and combustor geometry are coupled. 

Figure 3 shows a schematic of the flow where the inner separating 
streamline boundary may either enclose a recirculation zone or 
represent the inner wall of an annular combustor. Swirl is 
included in the domain outside of the recirculation zone. The 
various elements of the unified model are coupled together either 
through source terms or through the initial and boundary condi- 
tions. . _ . 

The features of the mixing analysis including assumptions 
are summarized below, (see also TABLE I) : 

The fuel droplets and other particulates form 
a dilute, continuum-like suspension in the 
multi-component gas phase carrier. 

. The volume occupied by the particulates is 
negligible. 

The particulates are classified in terms of 
size and composition and each class mixes at 
a rate determined by its ability to respond 
to the dynamics of the gas-phase motion. 
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. Turbulent transport is important only in 

the direction normal to the primary burner 
flow direction. 

Non-isotropy is accounted for by assigning 
appropriate transport coefficients to the 
axial and tangential components of the shear 
stress . 

The describing equations are parabolic and for steady 
turbulent flow are given as follows: 

Global Continuity : 
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where : 


N =0 for two-dimensional flow; N = 1 for axi symmetric flow 

\ • 

(w ). = production of the i^ 1 gas-phase specie due to 
^ 1 homogeneous gas -phase reactions; 

(w ) . = production of the i^ gas-phase specie from the 
P particle class due to evaporation or heterogeneous 

reactions; 

(Wp)^ = production of the jth particle class; 


Momentum Equations: 
x component 


du , du 
pu — + pv — 

3x dy 


_L d / N . . (rx)r v cm 3 .dp 

= — 7 —- ] y p (c ) a +r<r . B . — i - ~r" 

N 3y l v g L g . ] ]J 3y J dx 


(41) 


y component 
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; r = yw (circulation) 


(42) 


6 component 
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where for swirling flows only N = 1 has significance. 
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Energy Equations: 
mixture energy equ ation 
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where in the above describing equations the quantities 0^ , 6^ 
and CTj respectively define the ratio of the particle- to-gas phase 
eddy dif fusivi ties of mass, momentum and energy, i.e. , 
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and P , S c = gas phase Prandtl and Schmidt numbers , respectively . 
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Also, the gas phase mass fraction, a . , and the particle-phase 
mass fraction are defined such that 


a 

g 






where the subscript i covers each gas-phase species and the 
subscript j denotes each particle class. 


Note also that the energy equation. Equation (44) , may not, 
in general, be sufficient to define the thermal state of the 
system. Although dynamic equilibrium may be appropriate, the 
temperature associated with each particle class can be sensibly 
different from the gas-phase temperature depending on the process 
occurring on the particle scale. For example, for evaporating 
droplets the local saturation temperature is appropriate while 
if metallic particles are part of the multicomponent system and 
surface reactions are occurring then some other temperature is 
appropriate. Only in the case of inert particles (non- radiating ) 
is there similarity in particle drag and heating processes for 
which "near" dynamic equilibrium would imply near thermal equi- 
librium. In cases where the thermal equilibrium assumption does 
not apply. Equation (44) must be augmented by either the gas- 
phase or particulate -phase energy equation. 


In the present study the formulation for the mixing zone 
is limited to the thermal equilibrium assumption. Additional 
details of the treatment given to these equations can be found 
in Appendix A. 


Initial and Boundary Conditions 

The problem is set with the specification of the dependent 
property profiles at some initial x-station together with appro- 
priate boundary conditions invoked along the inner and outer 
"walls." 

These conditions are presented in TABLE I. There are, 
however, several implications associated with these boundary 
conditions as we have specified them: 

Inner Wall - When the inner surface bounds an imbedded, 
recirculation zone the state of flow along the surface is 
defined by the stirred reactor analysis. This state will 
generally be different from that in the surrounding directed 
flow. As a result diffusion across the dividing streamline ,y_ T7 , 

IW 
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occurs which provides the coupling mechanism between the inner 
and outer flow regimes. Of particular interest are the trans- 
port of species and energy. Since there is no net mass transfer 
across this boundary the respective species and energy fluxes 
out of the imbedded recirculation zone are given by: 


Species Flux, m. , 
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Energy Flux, e 
(rx) 


3a. 


e = 


. v g . rm _ Ml/1 _ s h 

Pr L 3y 3y i 3y 

J J 1 gas 


3/3 h 

P._E 

3y 


+ S /3 .6 . 


3h . 
_E1 


3 j 3y J 


+ E m.h . (T ) + E m ,h . (T .) 


1 1 g 


3 3 P3 


(47) 


Once the above fluxes have been determined over the entire 
dividing stream surface the total rate of transport for each 
specie and the energy may be found: 
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These specie and energy flow rates can now be (re) intro- 
duced into the stirred-reactor analysis providing an iterative 
method of bilaterally coupling the outer directed flow with the 
internal recirculating flow. Of course, a new average pressure 
level exerted upon the recirculation zone is also made available 
for updating the recirculation zone state. Finally, the contour 
of the recirculation zone must be specified and some data for 
basic swirling flows exits. Reference 16. 


Outer Wall - Lateral Injection - When the combustor flow 
encounters secondary injection through holes, the problem is 
reinitialized as described in TABLE I. The penetrations are 
determined by empirical correlations using the work of References 
17 and 18. Although several secondary injection analyses have 
pervaded the literature the aforementioned works were chosen 
arbitrarily, and the following discussion does not preclude the 
implementation of any other injection analyses. In fact, most 
of these empirical analyses treat the secondary injection problem 
in somewhat similar fashion. In particular, the interaction of a 
secondary jet and the subsequent dispersion in the primary 
(combustor field) fluid is assumed to be a two stage process - 
the penetration and mixing stages. In the penetration stage, the 
jet retains its identity while being accelerated and turned in 
the flow direction of the primary fluid. Here it is assumed 
that the jet emerging from the dilution hole acts as some "solid" 
body which is subsequently bent downstream by drag and distorted 
in cross section by pressure differences on the front and back 
faces and viscous shear. Accordingly, the jet trajectory using 
an analysis described in Reference 17 is obtained from 
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(50) 

where the quantity, q, defines the ratio of the dilution hole 
momentum flux to the primary stream momentum flux and the quan- 
tity, £x, defines the difference between the local axial 
distance and the axial location of the jet-dilution hole of 
diameter d.. Furthermore, the jet is assumed to deform from its 
otherwise Initial circular cross section into an ellipse having a 
5 to 1 ratio between its major to minor axis. Reference 18. 

The variation of the jet width or minor axis with distance is 
based upon subsonic experiments. Reference 17, and is given by 


h/dj = 2.25 + 0.22 (£x/dj) (51) 
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with the jet cross section area given by 

. A j (x)/A. (o). = ( 1 /5 ) (h/d ^ ) 2 (52 ) 

Thus, for given dilution-hole jet efflux conditions, the above 
three equations are sufficient to describe the penetration 
height, width and cross-sectional area once a (£x/dj) value 
has been chosen. Consideration of mass conservation and the 
fact that usually more than one dilution hole is located along 
the periphery of the combustor wall then allows us to deter- 
mine the geometric and physical properties of jet fluid effec- 
tively contained within an annulus having a radius as deter- 
mined by Equation (50) with a width given by 

*y/d. = (N./8) [(V&V'VaJ 

where N. defines the total number of dilution holes at the 
axial station in question. 

Once this annulus is defined and the jet fluid properties 
within it determine the mixing- analysis, or so-called second 
stage, is re-initiated. It is noted however that since the 
effective maximum penetration is assumed to be achieved quite 
quickly the distance that the jet traverses downstream of the 
port is considered to be small. Hence the second stage or 
re-initialization procedure occurs at the dilution hole axial 
station. 

Turbulent Transport Coefficients 

The method employed here to represent the effect of turbu- 
lence involves the specification of eddy transport coefficients. 
This implies a gradient mechanism of diffusion in direct analogy 
with molecular transport and this assumption is reflected in the 
form of the describing equations (38) - (44). 

Now, since the description of the . turbulence field is 
extremely complex a semi-empirical representation is required 
in expressing the eddy transport coefficients as functions of 
the state of flow. The presence of swirl further complicates 
the problem by introducing a tangential component of shear 
into the flow. 


V /a. r (53) 
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Although a complete description of the turbulent mechanism 
is not available there are some observations which are relevant 
to the present problem. First of all the effect of swirl 
enhances the axial component of shear and secondly, the shear 
distribution is anisotropic. References 19 and 20. 

These observations have been incorporated in the present 
turbulent transport property representations where the anisotropy 
is accounted for by specifying axial and tangential viscosities, 
respectively. 

Axial Viscosity, M ( rx ) — p C (tx) 

For small to moderate swirl where the axial velocity 
profiles are monotonic the work of Reference 21 has been found 
appropriate in representing available data. In particular, 
this model has the form: 

Mrx)= Ms) r; p (u ) (54) 

where k(s) is an empirically determined function of the swirl 
parameter, s, and is the half radius based upon the location 
of the mean velocity. Here, the data of Reference 20 has been 
used to define k(s). 


Values of k(s) found to represent non- reacting 
(isothermal) and flame data are as follows: 


. Isothermal runs 
k = 0.025 { s=0 } 

= 0.033 { s = . 1 } 
= 0.066 { s=.2 } 


. F lame 

k = 0.0225 
= 0.0225 
= 0*0450 


For large swirl, the axial velocity attains its 
maximum off the axis where the tangential component of the 
shear is zero. This suggests that the turbulent eddy scale 
on either side of the zero shear contour is different and 
that a two layer transport model is required. A limited 
study of such a representation was made using a model of the 
following type: 
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where 


y peak = y @ (pu) 

■ . max 


The results of this study are encouraging and comparisons 
with data are presented later. 


Tangential Viscosity, 




The previously cited work of Reference 19 has shown 
the anisotropic character of the shear distribution. This 
observation was based upon analysis of the measured axial and 
tangential velocity profiles indicating that the tangential 
viscosity is on the order of 10 to 20 percent of the axial 
viscosity. This conclusion is in general agreement with 
previous studies on vortex decay wherein a model of the follow 
ing type was suggested. Reference 22: 


M (re) = 
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(55) 


where r o is an initial value of the circulation (wy)init. 
and Vg is the gas phase molecular kinematic viscosity. 

Now, to complete the description of the turbulent 
transport properties, representations for the eddy diffusivity 
and eddy conductivity are required. For this purpose it is 
convenient to represent these coefficients in terms of the 
Schmidt No. and Prandtl No. respectively, viz: 


and. 
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Finally, the relationship between gas phase and 

particulate phase transport properties are expressed in terms 

of the parameters, a. , 9. and 6. , viz: 
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Although in the present study these parameters must be speci- 
fied some work on particulate dispersion in turbulent flows 
has been reported and is summarized in Reference 23. 

Now, the preceeding formulations have been applied 
to a range of configurations with particular attention given 
to the analysis of available data. 

4. MODEL APPLICATIONS 

To date the primary application of the model has been in 
terms of unit configurations involving both well-controlled 
laboratory combustors and actual systems. 

Pollutant emissions are most often measured at levels 
commonly considered as "trace" concentrations relative to the 
major products of combustion. Thus, their determination is a 
matter of great delicacy. Consequently, to establish a basic 
understanding of the mechanism (s) controlling the formation 
and disposition of pollutants has required well-controlled 
laboratory experiments. 

Now, during the course of the present study basic data 
was concurrently becoming available to us under a joint effort 
with ESSO Research and Engineering Company. In this regard, 
the work of Reference 12 should be noted. That study marked 
the first time that such basic systems as the H 2 /Air and Co/Air 
systems were studied experimentally and theoretically in terms 
of NOx emissions. A jet stirred reactor of the type developed 
by Longwell and Weiss, Reference 24, together with the above 
analysis were employed in that investigation. The results for 
the H 2 /Air and CO/Air systems are given in Figures 4 and 5 
where very good agreement between theory and experiment is 
shown for both systems. The significance of these results 
is two-fold:. First, the model can be used with confidence 
to accurately define the effects of variations in the con- 
trollable parameters. Insight gained from these two basic 
chemical systems can be applied to more complex fuels. 

Secondly, the H 2 /Air and CO/Air mechanisms are basic to any 
hydrocarbon/air oxidation mechanism. Thus > verification of 

these two systems provides a firm basis for postulating and 
developing schemes to represent the oxidation of hydrocarbons. 
Figure 6 shows a comparison of calculations using our quasi- 
global mechanism with data for the propane/air system. 

Included are equilibrium predictions and kinetic predictions 
also made assuming that free atomic nitrogen is present. 
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The latter is an artifice to establish the upper limit for 
the conversion of fuel "bound" nitrogen to N0 X . In addition, 
adiabatic as well as specified temperature results are presented , 
and the relative importance of the N+OHssNO+H reaction is also shown 

The adiabatic stirred reactor predictions are in good 
agreement with the fuel lean data but poor agreement is 
ohtained in the fuel rich regime. 

The calculations determined by using temperatures deduced 
from the experimental data yield low NO values throughout. 

The predictions are sensitive to small added amounts of 
free atomic nitrogen entering with the fuel. The conversion 
efficiency is high and can be traced to the fast N + O 2 = NO + 0 
reaction. 

The N + OH = NO + H reaction is most significant in the 
rich regime where [OH] > [o]» where [ ] is the concentration. 

Finally, the equilibrium results are as expected in the 
leaner regions but in the rich regime the measurements as well 
as the predictions of NO are higher than the corresponding 
equilibrium levels. This is apparently due to the dominance 
in this regime of free radical overshoots. 

In general, the above results are extremely encouraging 
and although more work is required in hydrocarbon kinetics 
modeling the results indicate that many of the observed trends 
are predicted with the quasi-global representation. Neverthe- 
less, our quasi-global mechanism does involve replacing many 
intermediate reactions and species high up in the oxidation 
chain with a subglobal partial oxidation step. The results 
presented so far have implied the importance of the non-equi- 
librium free radical overshoots on the generation of N0 X . 

In terms of the NO reactions considered in the present work 
0 , OH and N are crucial to NO formation mechanism. Of these 
species the combustion generated atomic oxygen appears to be 
the controlling one. For hydrocarbons, such superequilibrium 
values of atomic oxygen can account for an "early" formation 
of NO in the milli-second range that would not be predicted 
assuming the equilibrium level of atomic oxygen existed 
throughout. Figures 7 and 8 give an example of this for 
isothermal plug flow of octane and hydrogen assuming two models 
for the combustion process. In both models, the full system 
given in TABLE II was employed. However, the "coupled kinetics" 
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results were obtained by initiating the calculation with 
pure fuel and air whereas the "Equilibrium Fuel/0 2 - Non- 
equilibrium NOjj" results were obtained by initiating the 
calculation assuming equilibrium Fuel /02 wherein the N 2 
was assumed inert. Figure 7 shows that the atomic oxygen 
overshoot for the hydrogen/air occurs earlier than in the 
corresponding octane/air case. Of course, for pure air the 
growth of 0 is monotonic. The fact that hydrogen is more 
reactive than the hydrocarbon is the reason for "shift" in 
the overshoot and Figure 8 shows the consequences of this 
behavior. Thus, the early formation of NO is significant 
for the hydrocarbon but not so for hydrogen (or pure air) . 

The higher "final" pure air level of NO is due to the avail- 
able 0 2 that would otherwise be consumed by the fuel. 

Although this result does not cover a wide range of possible 
states that would be encountered in a practical system, it 
demonstrates, for a typical condition, that a coupled com- 
bustion/NO x mechanism is necessary to explain and predict 
certain observations. Moreover, the question that arises 
is whether any of the additional intermediates not included 
in the quasi-global mechanism play a significant role in the 
NO formation process where hydrocarbons are concerned. 

To gain some insight into this we directed our attention 
to the methane/air system. It seemed appropriate to start 
with this system since of the many hydrocarbons of interest 
this one appeared to be best understood in terms of its 
oxidation kinetics. The mechanisms used in this study are 
those given in TABLES III and IV. The results assuming 
adiabatic combustion are shown in Figure 9 where the data 
is from Reference 25. The jet stirred reactor previously 
cited was employed here. The Bowman-Seery mechanism pre- 
dicts the fuel lean NO emissions very well while poor agreement 
is obtained in the rich regime.. This behavior is remarkably 
similar to the results obtained for the propane/air system using 
our quasi-global mechanism. What is equally significant is 
the effect obtained when formaldehyde is included in our mech- 
anism. The comparison with the data is somewhat poorer than 
that obtained with the Bowman-Seery mechanism, although some 
improvement is observed in the fuel rich region. Although 
improved rate data for the formaldehyde chain may improve our 
predictions, these results support the contention that fast N 
producing reactions are missing as suggested by Fenimore, 
Reference 10. In any case, the common behavior between the 
quasi-global predictions and predictions made using "detailed" 
mechanisms suggests that the quasi-global concept itself should 
not be ruled out at this time. 
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There is an additional point regarding the effect of 
turbulence and the degree of "unmixedness " on the microscale 
in practical systems as well as within laboratory reactors. 
Studies on various aspects of this problem have indicated 
the potential importance of microscale inhomogeneities upon 
the combustion process. References 26, 27 and 28. In Refer- 
ence 28, Bowman, Pratt and Crowe introduce an empirical 
representation to account for incomplete mixing in their 
stirred reactor work. They showed parametrically that the 
temperature and NO emission levels are lower than would be 
predicted on the basis of assuming all fuel and oxidizer 
reacted within the geometrically defined combustion volume. 

This is, of course, the expected behavior when the fuel and 
oxidizer are premixed prior to injection into the reaction 
zone. In our work we have assumed a perfectly stirred con- 
figuration which for the Longwell reactor seems to be valid 
on the basis of the results obtained for the highly reactive 
H 2 /Air system. Figure 4. Now, in terms of the propane and 
methane results it is clear that if the model had assumed any 
level of inhomogeneity the predictions would have been displaced 
even further below the data. These results lend evidence to the 
observation that representing the kinetics rather than the 
effect of "unmixedness" is the most important consideration here. 

Now, all of the above comparisons with available data 
indicate that the models employed in our studies provide pre- 
dictions in complete agreement with the observed trends. 

These results have implied that insight into the effect of 
liquid fuel injection could be obtained and therefore might 
suggest some of the controlling parameters in droplet com- 
bustion. Accordingly, some calculations have been made for 
a situation where the fuel enters the recirculation zone in 
a partially vaporized state. The effect of residence time 
upon fuel-droplet consumption, recirculation-zone temperature, 
and NO concentration as a function of the droplet size at 
injection, are shown in Figures 10, 11 and 12. The complete 
set of operating conditions are given in TABLE V. Two values 
of the mass flow-to-volume ratio were selected to provide a 
range of residence times. The resulting residence times are 
about 2 msec and 17 msec. 

Figure 10 shows the importance of residence time upon 
droplet consumption. It is also interesting to note that the 
effect of droFwise combustion is small for both residence 
times meaning that the primary effect upon droplet consumption 


28 



is the heat conduction to the droplet. This effect is dominant 
for our conditions because of the relatively high temperatures 
in the surrounding gas phase due to the homogeneous gas-phase 
oxidation of the vapor-phase component of the fuel. These 
temperatures are shown in Figure 11 where it is interesting 
to note that at the larger residence time the amount of 
unconsumed liquid fuel is small enough for the entire range 
of: initial droplet diameters that dropsize has little effect 
upon the combustion temperature. Finally, Figure 12 shows 
the concentration of NO. The extreme sensitivity of NO to the 
temperature is clearly evident here. Of particular interest, 
however, is that at the larger residence time the NO levels 
are somewhat higher where the droplets are assumed to evaporate 
(with no dropwise burning) . This is a reversal in the trend 
observed at the shorter residence time. This is apparently 
associated with the fact that the difference in temperature 
between the two modes of droplet consumption is small. How- 
ever, the atomic oxygen levels are higher for the pure evap- 
oration mode of droplet consumption. Thus, under these con- 
ditions the important reaction: 

0 + N 2 = NO + N 

proceeds toward the right at a somewhat higher rate even 
though the temperature is slightly lower. This result 
further indicates the importance of the coupling between 
the combustion kinetics and the NOjj formation kinetics 
particularly where fuel droplets are involved. 

So far in our discussions emphasis has been on the emis- 
sions problem in terms of parameters controlling the fuel 
oxidation process. It is of interest, however, to examine 
other control measures of potential interest for NO^ reduction. 
In view of the preceding discussion involving liquid fuel 
injection it is appropriate to extend that analysis to include 
water injection here. 

Water injection is regarded as one possible control 
measure for NOx emissions. We have applied our multiphase 
stirred reactor model to the jet engine combustor data of 
Reference 29 and the comparison is shown in Figure 13. 

Excellent agreement is shown for the smaller size droplets 
and for the steam injection. The small difference between 
steam and the small droplet size predictions are in complete 
agreement with the experimental observations. The predictions 
also show that the effectiveness of liquid water injection 
decreases rapidly with increasing droplet size. Indirect 
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confirmation of this latter prediction has been. noted recently 
by Shaw, Reference 30, wherein he was able, upon becoming 
aware of the present results to explain the small reduction 
in NQx emissions on the basis of his dealing with large 
droplets. 

To further demonstrate the utility of this basic model 
element further comparisons with engine emissions data have 
been made. For this purpose, we have simply assumed that the 
primary zone is the source of NO and that the entire primary 
zone is a single perfectly stirred reactor. 

Figure 14 shows comparisons of predictions with measured 
NO x emissions from a variety of jet engines. The data as 
represented is a correlation as determined by Lipfert, Refer- 
ence 31. The predictions were made assuming an equivalence 
ratio of unity (<N1) and a constant residence time of about 
2 msec. Lipfert showed the importance of combustor inlet 
temperature with his correlation and the predictions are in 
complete agreement with this observation. The above comparisons 
of theory with data demonstrate that while certain questions 
on hydrocarbon combustion remain open much can be done with 
basic models to explain and quantify observations from practical 
system configurations. However, while this seems to apply to 
certain systems many practical burners depend upon simultaneous 
mixing and burning for effective operation. When reaction 
rates and mixing rates are of the same order of magnitude 
different modeling is required. In particular, most practical 
systems involve discrete injection of fuel and oxidizer which 
may involve imbedded recirculation zones and such systems are 
characterized by the existence of distinct diffusion flame 
regions. To analyze such flows the submodel which couples 
the reaction kinetics with the appropriate turbulent mixing 
processes is required. In these mixing regions it is assumed 
that convection dominates in the primary flow direction and 
that diffusive processes are most important in the direction 
normal to the primary flow direction as previously described. 

Initial applications of the mixing model analysis were 
made to non- reacting free swirling gas jet flows for direct 
comparisons of the results found in Reference (32). That experi- 
mental investigation involved a series of tests on free, non- 
reacting, turbulent jets with degrees of swirl covering the 
weak, moderate and strong ranges. The relevant data used for 
the comparisons found herein are given in Table VII. Measurements 



of mean axial and swirl velocities, static pressure and jet 
widths are reported at axial stations up to 15 orifice 
diameters. Figures 15, 16, 17 typify the agreement between 
the experiment and the analysis for the case where the 
swirl parameter, s, is equal to 0.6. 

Figure 15 depicts the decay in the maximum axial velocity 
for the so-called high and low swirl levels. Depending on the 
swirl level, the maximum axial velocity need not be located 
along the jet axis. In particular, for the large swirl case, 
s = 0.6, the axial velocity achieves a maximum off the jet 
axis for a distance of approximately 10 jet diameters. This 
can be seen by examining Figure 16 which shows the radial 
distribution of axial vddcity at several downs stations 

along with the theoretical predictions. The reason for such 
a behavior is that the presence of the swirl results in the 
setting up of radial and axial pressure gradients which, in 
turn, influence the flow field. In cases of strong swirl, 
the adverse pressure gradient along the jet cannot be overcome 
by the kinetic energy of the fluid flowing in the axial direc- 
tion and a recirculating flow is set up in the central portion 
of the jet between two stagnation points. This recirculation 
region extends downstream approximately 1.25 jet diameters. 
Reference 16. Corresponding to Figure 16, Figure 17 shows 
the radial distribution of the tangential, or swirl, velocities 
along with the theoretical predictions. The agreement is quite 
good within the inner region where the tangential velocity 
behaves like a solid body rotation? however, in the outer 
region the analysis somewhat under-predic ts the experiments. 

On the overall the agreement is quite good and although 
not shown fully similar comparative agreement has been indi- 
cated for the experiments involving the low and moderate swirl 
levels. 

Noteworthy also is that the theoretical results depicted 
have been obtained by utilizing a turbulent eddy transport that 
is both radially and axially dependent. x This model described 
in some detail previously was devised by examining the axial 
velocity profile histories for several degrees of swirl as 
well as the findings of Reference 19 where, through hot-wire 
measurements, it is shown that the nonisotropic character of 
the turbulent stress distributions increase with increasing 
swirl. The results of our theoretical efforts in attempting 
to establish closer agreement between theory and experiment 
have also indicated that the turbulent transport processes 
are strongly dependent upon the change in the character of 
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the flow as the level of swirl increases. For example, for 
low to moderate levels of swirl, i.e., for values of s < 0.6 
as adjudged by Beer and Chigier, Reference 16, the function 
k(s) found in Equation (54) was determined by deducing the 
value of eddy viscosity on the axis from the experimental 
profiles of velocity and pressure at a given axial station. 
Basically, the procedure was to solve Equation (41) for 
p(e v )g^rx) ( = n )at the axis where the derivatives of 
p and u required for this evaluation were obtained from the 
experimental profile information. 

For the non- reacting swirling flows, the results of these 
calculations are shown in TABLE VI. The same procedure was 
initially tried for the large swirl case (s = 0.6). However, 
it became readily apparent that one value of the function k(s) 
was not sufficient to provide a correspondingly comparative 
agreement indicative of the lower swirl results. In this 
connection it became clear that at least two mixing length 
scales and hence two levels of turbulent eddy viscosity are 
appropriate in the region of the flow where the velocity 
maxima occurs off the axis. In this case a preliminary two 
layer model, previously described, was utilized where the 
viscosity in the outer layer was predetermined to be related 
to the viscosity in the inner region which was deduced by 
the procedure described above. The values used for these 
comparisons areshown in TABLE VI as well. 

In addition, it was found that within the fully developed 
regions of the flow, the value of p, i.e., the ratio between 
^i( r ©) and , of between 0 to 0.3 gave essentially the 

same results. Whereas, in the initial regions of the flow 
where 1 the axial decay in tangential velocity and pressure is 
quite rapid, the agreement between the predictions and the 
experiment were more sensitive to the ratio. 

Accordingly, these basic studies which were performed 
primarily to establish the accuracy of the mixing model to 
represent the measured mean-flow properties in flows strictly 
governed by turbulent transport have clearly demonstrated 
that the non-isotropic character of the shear stress, alluded 
to previously, indeed had to be taken into account. Although 
good agreement was obtained with the available experiments, 
this required treating the turbulent transport coefficients 
with a certain degree of delicacy. 
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These observations notwithstanding, we have proceeded to 
make comparisons with swirling flame data. The experiments 
were performed using a swirl burner fueled by liquified petro- 
leum gas. Reference 33. The flame was unconfined and the 
measurements included velocity, pressure and temperature data 
along the flame. The results are shown in Figures 18 through 
25, and they include comparisons with data at two swirl levels 
where the assumptions regarding turbulent transport properties 
are given in TABLE VI. 

Figures 18 and 19 show the axial distributions of the 
streamwise and tangential components of velocity, respectively. 
The agreement is quite good at both swirl levels although a 
slightly slower decay in the predicted tangential velocity is 
observed in the downstream portion of the flames. Figure 19. 
Figure 20 shows a radial profile of the streamwise component 
of velocity and here also a somewhat slower decay is observed 
in the outer regions of the flow. Figure 21 shows the static 
pressure distributions along the flame axis. The agreement 
here is generally excellent- 

Perhaps of most interest are the temperature field and 
the NO x emissions from these flames. 

Figure 22 gives a comparison of the measurement and 
predicted radial temperature distributions for the large swirl 
level. Here the profile is just downstream of the flame 
closure point on the axis and the agreement is good. 

Figures 23 and 24 show the axial temperature distributions 
and again the agreeement is quite good although the measured 
peak temperature in the small swirl case. Figure 24, is some- 
what lower than predicted. 

Although pollutant emissions were not measured. Figure 25 
shows the emission index (E. I.) for NO obtained by integrating 
the NO flux across the flame at their respective closure points. 
The substantially longer residence time associated with the low 
swirl flame accounts for the higher NO emission level. 

Some comments are in order here with regard to both experi- 
mental procedure and the modeling of the problem. In the 
experiments, air was introduced within the swirl burner both 
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axially and tangentially with the degree of swirl changed by 
independently varying the axial and tangential flow rates. 

The fuel rich liquified petroleum gas/air mixture issued from 
the burner orifice and was subsequently ignited by means of a 
portable pilot flame. Accordingly, ignition was obtained by 
applying an "external" source which was removed upon establish- 
ing the flame. Thus, the calculations were performed by 
assuming a chemically frozen system from the burner exit to 
the experimentally determined flame stabilization point, i.e., 
x 0.085 meters. An equilibrium chemistry calculation was 
then performed in the neighborhood of where the stoichiometric 
ratio was unity to simulate the pilot ignition as it occurred 
in the experimental system. The products of combustion from 
this adjunct calculation provided an annulus of high- temperature 
gases which subsequently diffused into the inner fuel-rich 
regions and outer fuel-lean regions providing a simulation of 
the actual ignition source. 

It should be noted that for these reacting flow calcula- 
tions the function k(s) could not be determined from profile 
information as was described previously because there was 
insufficient profile information for these experiments. As such, 
the values of k(s) reported in Reference 19 were used both for 
the initial non-reacting region and subsequent reacting regions. 

In addition, some preliminary work was done in utilizing 
the mixing analysis to predict the flow field development within 
annular-type combusters. Figures 26 and 27 highlight some of 
the pertinent results obtained when the "centerbody option" as 
previously discussed (see also Appendix A) was exercised. 

Table VIII presents the initial data used with this option of 
the program and shows that the numerical experiment under con- 
sideration involves the subsequent mixing between air and a 
propane fuel jet issuing within an annulus at a speed of 
61 m/sec. The fuel is considered to be at ambient temperature 
while the air enters the annulus at a temperature of 1500° K. 

The initial pressure is assumed to be 10 atmospheres; and,, 
depending on the option employed the pressure either is held 
constant or an axial variation is allowed to develope as 
dictated by a prescribed area-contour. The duct outer radius 
is 0.328 m. (1 foot); the center-body radius being ^ ^ of the 
outer duct radius. The centerbody radius was held constant 
albeit the computer program can treat variable contour annuli. 

The fuel jet ring is located 0.254 m from the axis and has a 
ring width of 0.68(10) m, or approximately 0.3% of the total 
width of the annulus. 
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The relatively small width of the fuel-jet ring brings 
out another interesting aspect of the flexibility of the 
computer program and that is with regard to mesh size required 
for intial profile definition- For example, if one wishes to 
input into the program, data across the entire annulus while 
still maintaining good profile resolution between the airstream 
and fuel stream an inordinate amount of input information, due 
to the small grid spacing needed, would be required. However, 
because the flow field bounding the jet initially is uniform 
all that is required is detailed information in the regions 
just surrounding the jet. Accordingly, as mixing proceeds, 
grid points are added from above and below the fuel- jet center- 
line as required to bring into proper perspective the necessary 
definition of the potential flow regions. Thus, as indicated 
in TABLE VIII, information at only eight grid points are supplied 
initially; four points to define the fuel jet and four points 
(two on each side of the fuel jet) to define the surrounding air 
stream. Further details pertaining to this option can be found 
in Appendix A. 

Figure 26 shows the downstream development of the tempera- 
ture within the annulus due to mixing of the high- temperature 
air stream with the cold fuel jet. These results were based 
upon imposing a constant wall pressure (10 atm) along the duct 
and allowing the outer wall to adjust to this constant value of 
pressure . 

The numerical results show that the outer- wall contour did 
not vary markedly within the 1 radii of duct length considered. 
The figure also shows that the mixing process has reached the 
inner, center- body wall before spreading across the entire duct 
width. This asymmetry is due primarily to the closer proximity 
of the fuel jet to the inner wall than the outer wall and because 
the eddy- viscosi ty value is assumed constant throughout the duct. 

As indicated previously, the other option investigated was 
one in which the outer-wall contour is prescribed a priori ' 

(equal to 0.328 m. , and constant) and hence the pressure field 
within the duct is adjusted at each diffusion step to reflect 
this constraint as well as fulfilling the requirement of mass 
conservation. Figure 27 shows the tempera ture- prof ile develop- 
ment under these conditions and when compared to Figure 26 shows 
the effect of the adverse pressure gradient that arises due to 
viscous mixing. In the previous case, no pressure gradient 
is set up because of the imposed constraint of constant wall 
pressure. 
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Overall, these results were encouraging from the stand- 
point of having established a model capable of lending insight 
into unit -problem behavior. 

In general all of. the above studies have demonstrated that 
much can be done with basic models to aid not only in understanding 
the parameters controlling flame -genera ted pollutants but also in 
defining areas in kinetics and turbulent mixing needing additional . 
experimental and theoretical work. 



5. 


SUMMARY OF RESULTS 


The study performed here has resulted in the develop- 
ment of some of the basic elements required in turbojet 
combustor modeling. Application of subelements of the unified 
model has provided needed insight into the NC>x formation pro- 
cess while also establishing the validity of many of the 
included mechanisms. In particular, the following constitute 
the major observations of the study: 

1. The kinetics mechanisms for the I^/Air and CO/Air 
systems appear adequate for combustion and NO emissions 
predictions for temperatures at least in the 1500 to 2200° K 
range . 

2. The quasi-global mechanism underpredicts the 
observed NO emission levels in the fuel rich region. 
Preliminary studies with detailed methane mechanisms show 
the same behavior suggesting that the quasi-global concept, 
is not at fault but rather certain fast N producing species 
and reactions are required. 

3. Preliminary studies with liquid fuel injection and 
comparisons with data for the effect of water injection on 
NO emissions show the strong influence of droplet size. 
Excellent agreement with water injection data was obtained. 

4. All studies made here show NO 2 levels to be from 
two-to-three orders of magnitude smaller than NO. 

5. The stirred reactor concept with appropriate kinetic 
mechanisms explains much of the observed data obtained from 
both laboratory reactors and practical engine configurations. 

6. The mixing model predicts the observations made on 
non-reacting and reacting swirling flows but better agreement 
between theory and experiment requires treating the transport 
properties with great delicacy. More work in required here. 

7. The unified model has been completed but its appli- 
cation and further development are required. 

In general, the study not only has provided much needed 
insight into the NOx formation and disposition mechanisms but 
has also delineated areas in the kinetics and fluid mechanics 
needing additional work. 
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6. APPENDICES 

APPENDIX A - MIXING EQUATIONS IN 
VON-MISES COORDINATES 


The starting point for the coupled mixing and kinetics 
model is the boundary-layer form of the conservation equations 
for global mass, momentum and energy and species diffusion, 
i.e.. Equations (38) through (44). A solution of this system 
provides the details of the flow field including axial and 
tangential velocity, pressure, temperature and species fields. 

Now, by introducing a stream function such that 

M N 

0 . d0/3y = puy 

and (A-l ) 

M N 

0 30/dx = - pvy 

global continuity is automatically satisfied and the resulting 
coordinate transformation from the physical (x,y) system to a 
(x , 0 ) system yields a set of differential equations which are 
more amenable to finite-difference solution techniques. Appli- 
cation of this generalized stream function yields the following 
set of parabolic, partial differential equations : 
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Mixture Energy Equation 
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Once initial and boundary conditions are specified, and 
assuming thermal equilibrium for the particulate phase, the 
above set of equations, along with the equation of state, i.e.. 



p/ { RT ( L OL / w i + S )3 j / W /) } 


(A-8) 


together with thermochemical data for the specific enthalpies 
in the form 



(T) 


for the gas-phase species and 


h . = h . (T) 

3 3 

for the particular-phase species and expressions describing 
the volumetric' productlo n rates of each specie properly define 
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the problem. In the above system of equations the quantities, 
J, reflect the interaction between the various phases within 
this multiphase system. Accordingly the J's are defined as 
follows : 
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(A— 10) 

(A— 11 ) 
(A— 12 ) 
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(A-14 ) 
(A-l 5 ) 


Also, for simplicity ( the parameter 6. appearing in Equation (A-7) 
is defined by J . 

e. = -Le0 - (h /h.f: /(Le - 1) 

. * j y j 

An interesting feature of the analysis written in this form is 
that it suggests an appropriate modification of the gas-phase 
transport coefficients making it applicable to "near" dynamic 
equilibrium multiphase flows. For example, Equatiors (A-4 ) 
and (A— 11 ) • show 
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is the gas-phase turbulent viscosity. Thus, assuming that 
the functional form of the gas-phase transport coefficient, 
i.e., Equation (54) is unchanged due to the presence of 
the particles one can apply an appropriate model for fi (rx) 
and arrive at an effective global model for |i T according to 
the above equation. 


Likewise, Equations (A-9) and (A-10) can be thought of 
as being used to define the "effective" gas phase and par- 
ticulate phase Schmidt numbers in that one could define 

< 1/Sc) g,eff = (V s c)g ri ‘ . 

(1/Sc) p,eff '= e /< Sc >g . 

These equations can be readily put into a finite-difference 
form where "backward-difference" scheme is employed for the 
axial derivatives and where central-difference scheme is 
used for the derivatives in the radial direction. But, 
before this is discussed, the various boundary and initial 
conditions which one can implement with the above analyses 
will be documented. 

As mentioned in the body of this report several types 
of jet configurations can be analyzed with the above mentioned 
scheme. These included: 

free jets 

ducted jets 

ducted annular jets 

For the latter problems the annulus could be either a 
physical center body, which is relevant to combustors or 
the "center-body" might be the outer boundary of a recir- 
culation zone. Although, from a geometric standpoint both 
of these submodels are similar, the boundary conditions along 
the "center-body" for each must be treated differently as 
will be discussed later on. 
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For the initial conditions for each of these models one 
must specify the distributions of velocities, species, pres- 
sure and temperature, i.e., V 

x = o : 

u = u (o,0 ) 

r = r(o,<M 

p = p (o , 0 ) * 

ot. = ot (o,0) (A- 18 ) 

#i = ^i (°,^ ) 

T = T(o ,0)** 

In addition to this standard set of profile information 
one requires an initial distribution of 

dy/dx = dy/dx (o ,0) 

dp/dx = dp/dx (o ,0) (A- 1 9 ) 

0 j , 6 j ,ct j = 0j ,6 j ,o\(o,0) 

u = M (o , 0 ) 

where the latter two forms of profile information are assumed 
to be unchanged as the mixing proceeds. 

As mentioned previously the boundary conditions for the various 
jet configurations must be handled separately depending whether 
an axis of symmetry or a centerbody exists. Thus the boundary 

* Note: if the initial distribution of tangential velocity 

or r(r = wy) is specified then p(o ,0) can be 
determined by quadrature (see Equation (A-5 ) ) 

** Note: Although one of the dependent variables is the 

total enthalpy, H, the idea of requiring an initial 
temperature distribution instead is totally con- 
sistent since H(o,0) can be determined from this 
and the other profile information. 
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conditions for the free-jet mixing properties are in effect: 
Free Jet Boundary Conditions 

0 = 0: ! bu/bll) = aOj/90' = a£_./50 = = r = o 

. ... CD 
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The governing equations, i.e. , Equations (A-2) - (A-7) 
must be treated differently along an axis of symmetry. 

For as the axis is approached 

. 0-0 

and 

'd/dlp - 0 


Accordingly, by applying L'Hospital's rule, the governing 
equations applicable along an axis of symmetry are 
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For ducted jets with no "center-body" the conditions of axial 
symmetry, i.e.. Equation (A-20) still prevail but the free- 
stream conditions that make up the outer boundary conditions, 
i.e.. Equation (A-21) for the free jet problem must be replaced 
to reflect the fact that the outer boundary is now the wall of 
the duct. 

Ducted Jet (no centerbody; outer contour prescribed) 


0 = 0: same as Equations (A-20) 


(A-22 ) 


0 = 0 


wall : 


ciOj/90 = 9£_./d0 = 0 
3H/B0 = 0 


impermeable wall 


T = T 


wall 


adiabatic wall 


— — prescribed wall 
temperature 


(A-23) 
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*/*-- (* M /Pw N ) wall <P“ 2 /t« g lrx) ) mean : 

a<r/yVi>* = - (* M /puy N+1 ) ua i 1 (P u3 /M g <rx) ) mean (T/wJ, all <c£/2) 


The latter two conditions express the fact that in reality 
a boundary layer is developing along the wall. Thus, in 
effect two length scales define the mixing region, namely 
the boundary layer thickness along the wall and the radius 
of the duct. Since, the boundary layer thickness is very 
small compared to the duct radius, obtaining its detailed 
flow structure near the wall would require a highly refined 
numerical mesh relative to the mesh needed for good resolution 
of the bulk flow field. It is assumed here that the details 
of the wall boundary layer are not a dominant influence on 
the development of the bulk flow field and in this connection 
the gross effects of the wall boundary layer are included by 
using these two equations without resorting to unwarranted 
detail. Accordingly, the boundary conditions on the velocities, 
u = 0, w = 0 are replaced by .relations between wall shear and 
velocity gradients in terms of the skin friction coefficients 
whose values are considered as input into the numerical programs. 
Thus, these equations can be derived as follows: 

T <r*> /M <rx> = _ ( 3u/8y , wall = <' ^,» <pu* >_„/* 


■: T (r 8 ) /u (re) = - [ s (r/y s )/s y ] 
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wall 
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* Note r and y both reflect the radial coordinate in the case 
where N = 1 . 
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For the ducted annular jet, the centerbody is considered as 
an impermeable wall. As such the outer boundary conditions 
associated with the previous problem also apply here for 
the inner wall. - 



sr/yVa* - W M /puy H+1 ) IW <puVu g (Kt) ) mean <?/>) IW < c f/2) IW 

There is a difference however between these two problems 
and that is for the centerbody problem it is possible for 



which can yield a non- regular behavior along the centerbody 
if the exponent, M, in Equation (A-l) is other than zero as 
evidenced by examining the inner wall boundary conditions. 

For example, if M = 1 as is usual for treating free-jet flows 
Equations (A-2 6) for u and r would show that du = ^(T/y 3 )/ 0 

even if c f/2 / 0! This difficulty is circumvented by imposing 
M = 0 and therefore shows the utility of the generalized stream 
function as defined by Equation (A-l). 


47 



Now, if the so-called centerbody is considered as the 
outer bounds of the recirculation zone then the "inner wall" 
boundary conditions must reflect the fact that across the 
boundary there can exist a diffusive transfer of mass and 
momentum as well as energy. In this regard then what is 
specified initially along the "dividing streamline" are the 
axial distributions of the dependent variables as determined 
by the well-stirred reactor analysis. Thus, as a first iterate., 
interaction across the dividing streamline will be unidirectional 
and the kinetics associated with the recirculation zone will 
affect the processes occurring within the jet mixing region 
whereas the reverse situation is, during this first iterative 
"pass," not considered. In essence then, the recirculation 
zone is treated as a passive source (or sink) of thermal 
energy and species . Once this first "pass" is made, evaluation 
of the fluxes of mass, momentum and energy across this dividing 
streamline, as discussed in the body of this report, will 
define a new recirculation zone state for treatment with the 
well-stirred reactor analysis and thence the mixing analysis 
can be repeated. Accordingly the boundary conditions for this 
sub-problem are: 

Ducted Annular Jet (Recirculation Zone Contour Prescribed ) 


0 = 0 : 
0 = 0 


OW: 


Same as Equation (A-20) 
Same as Equation (A-23) 
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Essentially this completes the analyses required for 
all the options and submodels associated with the mixing 
model. However, as TABLE I indicates, the ducted jet. sub- 
model can also treat problems where instead of providing 
a prescribed outer wall contour one can impose a prescribed 
pressure. In this connection, the unknown contour is 
deduced from the calculations by imposing conservation of 
mass at each axial station. 

In addition, the equations, as set forth here, were 
grouped for ease in the discussion and in formulating their 
finite-difference counterparts. However, the actual computer 
program logic need not reflect blocks of calculations 
commensurate with the internal groupings within the equations 
as are shown here. 



APPENDIX B - FINITE -DIFFERENCE FORM OF GOVERNING EQUATIONS 


The governing partial differential equations describing 
the conservation of mass, momentum, and energy within the flow 
field and along the axis can be readily put in an explicit 
finite-difference form. A backward-difference scheme is employ- 
ed for the axial derivatives and a central-difference scheme for 
the radial derivatives. 


Consider the flow field 
divided into a grid in 0, 
x coordinates . Then the 
derivatives of an independent 
variable, say F, will be evalu- 
ated by 
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where 


a , i = (a + a )/2 

n,m+^ n,m n,m+l 


(B-4) 


The M+l-ordered derivatives required for the equations 
along the axis are put into the form 
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Accordingly, using the formalism as described by the above 
equations, we have 
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A A*) 2 


$ 


•IT (A) (— ) b (l-Jr. )] < 0 ^> } ( B -7 ! 

IL Pr g • D . Jn,m n,mJ 


( Ai) , . = Ax 

n+l,m 


r ( dY) (J£-) _ (^E) A) } + 

l dx n+l,m y 3 u n,m dx n+l,m pu n,m ■> 0 M (&0) : 


{ (A) J (x) <u> 

L L m J n , m n , 


m J 


(B-S 


Ur) n+l,m* ITT ^ { Ca A J M le> 


^ L~(rx) ; 


0 Ai/j 


n,m |n,m M ( ^ |rM 


0 A0 


- 1 n,m 
( B- 9 ) 
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AJ_ 


(AH) 


f A(L -1) : 


4* { f — <h> !* + * { r . ^ 

n+l,m M . 2 1 ( 1 - P.. J _ _ n,m J ( 4 j P_ I n,m l J 


h'.<or 1 > f- 


0 t*l)‘ 


r n,m 


i r 


*{ ?[ 


A ( L -1)_ 


P 

D ~ r 


0 h . < 0 .> 
n,m j 3 D 


n,m ■/ 


+ * jTaU • (JJ.-.1) - r 0 6 .M./u)”. <q a /2> } 

lL g g j J 3 J J n,m J 


o rr ^(q, (vi-D - s j?.o-.Myu)i <w a /2> 

XL g • 3 3 3 g j n,m n,m J [ 


Ax ( 


T ( g ). 

h o' 

d g q 1 


(p) 


0 


^ * A — j- +$ j A £ 0.h.0. j 

• v 1 a J n,m I s 3 j j Jn,i 


■VI ’ 1 s 


+ * Jm (8) m } m . 


m 


(B-10) 


The radial momentum, i.e.. 


M 


?p /*$ = >/) t^/uy 


( B— 11 ) 


can be obtained by simple quadrature and the quantity A in the 
above equations is defined as 


A . 2N ( rx ) / , M 

A = y u pu/w, 

g 


( B-12 ) 


Application of the finite-difference forms, i.e.. Equations B1 
through B4 are also applied to the requisite boundary conditions 
that have been discussed in Appendix A. 
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The axial step size. Ax, must be kept small to ensure sta- 
bility. This is done by setting Ax equal to the smallest of the 
triplet (Ax^, Ax 2 » Ax 3 ) where the following criteria is used to 
establish their values: 


(1) For 


m=2 


V 55 V 

m=Me 


Ax, = 


,N /Af xM+l 

ife iMl 


(— ) + (—) 

S n,m+h Sc n,m-^ 

c 


( B-13 ) 


( 2 ) 


For & = >J) 


m=l 


Ax 2 = 


IM1 


M+l 


(1-N) 


3 I 2 


i' ■(n+ik , ,n 


{0 U) (£-) . 

n,m Sc n, 1 


( B-14 ) 


(3) Ax_ = minimum (Ay) (B-15) 

3 m 


Then Ax = min ( Ax^, Ax , Ax^) and accordingly Ax is directly 
related to A:0 (an input) and inversely related to the local 
value of eddy viscosity. 
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APPENDIX C - NOMENCLATURE 


A factor which differentiates between diffusion-- 

controlled combustion (A=l) and droplet evaporation 
(A=0), Equation (24) 

A. cross sectional area of lateral jet; dilution hole 

-i 

J analysis. Equation (51), m s 

B transfer number in droplet combustion analysis. 

Equation (24) 


skin-friction coefficient 

C . v mole fractions 

ID 

specific heat, J/kg/°K 

d droplet diameter, m 

jet diameter, m 
dilution hole diameter, m 

activation energy, k cal/mole/°K 

total rate of energy flux across dividing streamline. 
Equation , (49) , j/k9 

energy flux across dividing streamline/unit surface 
area. Equation (47 ) f j/kg/m 3 

heat of combustion. Equation (24 ) t j/kg 


total enthalpy of mixture, j/kg 

enthalpy of gas-phase mixture in multi-component system, j/kg 
static enthalpy of specie k , J/kg 

mixture enthalpy or width of lateral jet. Equation (51), J/kg 
defined in Equations (A-9) through (A-15) 
k thermal conductivity, J/ m sec °K 


k,_ , k, forward and backward reaction rates for reaction p 

fp bp 
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k 

°: f ,P 

k (S );' 


Le 

«T 

m 

N 

N . 
D 

N 

P 

Pr 

Q 

Q 

q 

q 

R 


s 

Sc 

T 

t 


reaction rate constant for reaction p 

empirically de termined cons tants used in eddy^ viscosity 
models of mixing subprogram. Equation (54) 

Lewis number 

total rate of mass flux across dividing streamline, kg/sec 

mass flow rate within stirred reactor, kg /sec 

geometric parameter, N=0 implies 2-D flow, N=1 implies 
axisymmetric ..f low , or number of gas-phase species 

total number of dilution holes at an axial station 
or numbe r dens i ty 

stoichiometric oxygen/fuel ratio 

pressure , N/m 2 

Prandtl number 

net rate of external heat addition 
heat of vaporization 

ratio of dilution hole momentum fl\nc to primary stream 
momentum flux. Equation (50) 

component of velocity in x,r plane, m/sec 
gas constant, j/°k/fool 

jet "half radius " based upon location of the mean 
velocity within mixing region, m. 

sWirl factor, = angular momentum flux/axial momentum flux 
times nozzle radius. 

Schmidt number 

temperature, °K 

time, sec “• 
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u,v,w components of the velocity vector in the axial (x) , 

radial (y or r) , and tangential (8) directions, m/sec 



V. 

D 


W_ 


W. 


(W g’i 




(w ) F 

P 3 


(y i> g 


IW 


volume of stirred reactor, m 3 

.average velocity of lateral jet. Equation (53), m/sec 

molecular weight of specie k, kg/kg mQ ^ e 

volumetric rate of production of specie k, kg m 3 /sec 

production of the i tn gas-phase specie due to 
homogeneous gas-phase reactions, kg/m 3 /sec 

production of the i^ 1 gas-phase specie from the j th 
particle class due to evaporation or heterogeneous 
reactions, kg/m 3 /sec 

production of the jth particle class, kg/m 3 /sec 

mass function of species i in the gas-phase subsystem 

radius of dividing streamline of innerwall of 
annular combustor 


y radius of outerwall of combustor, m 

(Ay ) "Effective" ring width due to dilution hole efflux, m 


a mass fraction of gaseous specie k 

k 

a mass fraction of all gas-phase species in multi- 

^ component system 




/ 3 . 


mass fraction of particulate species k 

mass fraction of all solid (or liquid) phase species 
in multicomponent system 


r circulation, T = wy , m 3 /sec 

6 droplet or bulk density , kg/m 3 

€ D #€ v' C T turbulent gas-phase exchange coefficients of mass, 
momentum and energy, respectively 

8., 6 .,cr. ratio of partic le-to-gas phase eddy diffusitives 
■* ^ of mass, momentum and energy 
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X 


burning or evaporation rate constraint. Equation (23) 




eddy viscosity , N- sec/m 3 






// 

V ip 


P 

T 

T 

$ 

V 


ratio of eddy viscosities in the rx-plane to the 
r 0-plane 

molecular kinematic viscosity, N-sec/m 3 

stoichiometric coefficients of reactants i' in reaction 

stoichiometric coefficients of reaction products i in 
reaction p 

density, kj /m 3 

residence time within stirred reactor, sec 


ratio of shear stress in (rx) to (r0) planes 
equivalence ratio = fuel^air ratio/<p 
stoichiometric fuel-air ratio 




stream function 

Subscripts 

[ 

] 

molar concentration, moles/cc 

( 


implies liquid state 

( 

> 

implies gaseous state 

( 

} o 

conditions along axis 

( 

) 

ow 

outer wall 

( 

) iw 

inner wall or dividing streamline 

( 

} p 

implies particulate state 

( 

)„ 

free stream conditions 

( 


maximum value 

Superscripts 


( -) 1 

( )° 

( ) B ' P * 
, >rx 


implies inflow state for Recirculation Zone Model 

implies outflow state for Recirculation Zone Model 

boiling point 
•rx-component of ( ) 
rf-component of ( ) 
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TABLE' l. ~ COMPUTATIONS IN MIXING ZONES 


I 

VARIABLE MIXING . , 

ZONES 

11 

SPECIFIED INITIAL 
CONDITIONS 

III 

BOUNDARY CONDITIONS ALONG 
IMPERMEABLE PART OP OUTER 

IV 

BOUNDARY CONDITIONS ALONft INKER WALL* 

A\ l a i Velocity, u(x,y) 

“ * ufo.y) tor y lw • y •: 

S ye'-* 1 1‘ , 'ati sKin friction 
of shear: 

c.(x) , » 

f rx or r r* * 

rMur.*.' ar 
Flo- 

Q^ibeddod 

Recirculation 

Zone 

Same as outer w.*il 

Specify i: JW » u^U) hated upon 

mean flow tJ:r;*>*.gj velocity j.r# 
i ecirc’ i* t L:n . no 

Swirl Velocity, w(x,y) 

w * w(o.yl for y JW ' y ■' y w 

Specify «**aii akin friction 
or shear: 

c.(x) . . 

f re or t xb m 

©0 

Same outer «»*aM _ r ^ 

rperifv - Iw -v lw «•>./„, 

Tenpe * 4 tore, T(x.y) 

T » T(o,y) for y w *; y < y w 

Specify adiabatic or wall 
temperature: 

• 0 « v*> 



Same it outer wa > 1 

@ 

Specify H - H reel rcu* at ion 

rone 

Maas Fraction, O, ^(x.y) 

Number density. N^ (x.yj 
Droplet (Particle) air* d .(x,y) 

“i.J * °i.i {o - Y) 

tor y, w < y < y w 

"j • V°' V,,or y w < y < y^ 

<J j - dj (o.y)for y IW < y < 

Impe rmeabi 1 i ty ; 

^4 • 

v* ■ ° 

© 

Same as outer wall 

© 

Specify recircu- 

lation lone 

Pressure , p (x,y) 

p - p (o.y) for y llf < y < y w 

Specify y QW (x) and p^U) 

is obtained as part of the 
solution gjr obtain y^tx) 

if P w (x) la specified 
(9ee Wall Contour below) 

© 

Determined as part of the 
solution 

© 

Determined as part of the 
solution 

Mall Contour Outer 

Inner 

<« 

if y ow* peclfi * <s * hen 

la obtained as part of the 
solution ar obtain Y w l*) 
If ia specified 

© 

Sped fled i 

y iw ** y iw (x> 

»» -yp. <°> 

© 

Specified: 

y j ¥ * yiK 1 "' 




V 

AIR INJECTION ALONG THE WALL 


Variable 

Slots 

Holes 

v r 

T i 

”i 

A 

When the combustor flow reaches' a slot, 
the problem ia reinitialised. 

The "new" Initial conditions involve 
the extension of the prof lee at x x , 
as given in Column II, across the slot 
where the slot width ia known by the 
•peel flea tion of * y^ix). 

When the combustor flow reaches 
a hole, the penetrations are 
determined and the’ flow properties 
at x - S£ are specified as shown 
In Column 11. 

For both slot* and holes, the state 
of the air flow (f^ , Tj * Vj) must 

be specified on the upstream side of 
tii« openings. 7T»en, given an estimate 
for discharge coefficients, the state 
of flow entering the combustor is 
determined. The distribution of slots 
and holes are input infomutirn. 


for full flow aynnctry conditioxi «r« applied on I ha axis 




TABLE II 



EXTENDED C- 

■H-O CHEMICAL KINETIC REACTION MECHANISM 




X £ » AT* exp(- E /RT) 

FORWARD 



reaction 

A 

b 

E/R 

1. 

C H 0 - j^H +n CO 

n m 2 2 2 2 

5.52x10® _«« 

- . C C 

.825 C H 0. 

p n m 2 

1 

12.4 x 10 3 

2. 

00 + OH » H + CO^ 

5.6 x 10 U 

0 

.543 x li 

3. 

00 + o 2 > co 2 + 0 

, ,„12 
3 x 10 

0 

25.0 x 10 3 

4. 

00 + 0 * M = 00 2 + M 

19 

1.8 x 10 

-1 

2 x 10 3 

5. 

H2 ♦ 0 2 = OH + OH 

1.7 x 10 13 

0 

24.7 x 10 3 

6. 

OH + H, » H, 0 + H 
2 2 

2.19 x 10 13 

0 

2.59 x 10 3 

7. 

OH + OH * 0 + H 2 0 

5.75 x 10 12 

0 

.393 x 10 

8. 

0 + H 2 = H + OH 

1.74 x 10 13 

0 

4.75 x 10 3 

9. 

H O «• O + OH 
2 

2.24 x 10 14 

0 

8.45 x 10 3 

10. 

M + 0+ H=>0H+M 

1 in 16 

1 x 10 

0 

0 

11. 

M + O + 0 = 0 2 + M 

9.38 x 10 14 

0 

0 

12. 

h + h + h»h 2 +m 

c ,« 15 

5 x 10 

0 

0 

13. 

M + H + OH = H 2 0 + M 

,„17 
1 x 10 

0 

0 

14. 

O + N 2 - N + NO 

1.36 x 10 14 

0 

3.775 x 10* 

IS. 

N 2 + °2 " N + N °2 

2.7 x 10 14 

-1.0 

6.06 x 10 4 

16. 

Nj + 0 2 . NO + NO 

9.1 x 10 24 

-2.5 

6.46 x 10 4 

17. 

NO + NO * N + NO 

2 

1.0 x 10 10 

0 

4.43 X 10 4 

18. 

NO * O = 0 ; + N 

1.55. x 10 9 

1.0 

1.945 x 10^ 

19. 

M+NO = 0 + N + M 

2.27 x 10 17 

-0.5 

7.49 x 10 4 

20. 

M + N0 2 = 0 + NO + M 

1.1 x 10 16 

0 

3.30 x 10 4 

21. 

M + N0 2 » 0 2 + N + M 

6.0 x 10 14 

-1.5 

5.26 X 10 4 

22. 

NO + O, = NO, + O 
2 2 

, ,„12 
1 x 10 

0 

2.29 X 10 4 

23. 

N + OH = NO + H 

4 x 10 13 

0 

0 

24. 

H + trc> 2 s> NO + OH 

3 x 10 13 

0 

0 

25. 

00. + N o CO + NO 
2 

2 x 10 11 

-1/2 

4 x 10 3 

26. 

CO 4- N0 2 = C0 2 + NO 

2 x 10 U 

-1/2 

2.5 x 10 3 


Reverse reaction rate, , is obtained from ^nd 
the equilibrium constant, K . 



TABLE III 




CH /0/N 

System 

(Reference 14 , 

Bowman/Seery) 





k f 

= AT^ exp ( - E /RT ) 




REACTION 



FORM ART) 






A 

b 

e /r 

1. 

CH 

4 

+ M = CB 3 +H+M 


, 17 
2 x 10 

0 

44.5 x 10 3 

2. 

CH 

4 

+ OH = CH 3 + H 2 0 


2.8 x i0 13 

0 

2.5 x 10 3 

3. 

<*4 

+ O = CH 3 + OH 


2 x 10 

0 

4.64 x 10 3 

4. 

<*4 

+ H = CH + H 
3 2 


6.9 x 10 13 

0 

5.95 x 10 3 

5. 

<*3 

+ ° 2 - h 2 ° + CHO 


„ ,.10 
2 x 10 

0 

0 

6. 

CH 

3 

+ O = CHO + H 

2 


1 x 10 14 

0 

0 

7. 

CHO 

+ OH = CO + H 2 0 


1 x 10 14 

0 

0 

8. 

CHO 

+ M = H + C0+M 


2 x 10 12 

1/2 

14.4 x 10 3 

9. 

0 + 

N 2 = NO + N 


1.4 x 10 14 

0 

3.79 x 10 3 

10. 

N + 

0 2 = NO + O 


6.4 x 10 9 

1 

3.14 x 10 3 

11. 

N + 

OH = NO + H 


4 x 10 

0 

0 


64 



TABLE IV 



GASL 

CH 4 /0/N System 




k f- 

AT b exp(-E/RT) 




REACTION 


FORWARD 




A 

b 

e /r 

1. 

CH + M = CH + H 
4 3 

2 x 10 

0 

44.5 X 10 3 

2. 

CH 4 + OH = CH 3 + H 2 0 

3.5 x 10 14 

0 

4.5 x 10 3 

3. 

CH + O = CH + OH 
4 3 

2 x 10 13 

0 

3.45 X 10 3 

4. 

CH, + H = CH +H 
4 3 2 

2 x 10 14 

0 

5.95 x 10 3 

5. 

CH 3 + 0 = HCHO + H 

1.9 x 10 13 

0 

0 

6.* 

CH, + 0, = CHO + HO 
3 2 2 

2 x 10 10 

0 

0 

7. 

CH 3 + ° 2 = HCHO + OH 

1 x 10 

0 

0.75 X 10 3 

8.* 

CH 3 + O = CHO + H 2 

1 x 10 

0 

0 

9. 

HCHO + OH = CHO + H 2 0 

3 x 10 

0 

0 

10. 

HCHO + H = CHO + H 2 

1.7 x 10 13 

0 

1.5 x 10 3 

11. 

HCHO + CH = CHO + CH 
3 4 

2.5 x 10 10 

1/2 

2.65 x 10 3 

12. 

HCHO + 0 = CHO + OH 

3 x 10 13 

0 

0 

13. 

HCHO +■ 0, = CO, + H,0 
2 2 2 

7.3 X 10 10 

1/2 

0 

14. 

cho + o 2 = oo 2 + OH 

7.4 x 10 11 

1/2 

0 

15. 

CHO + O = <X> 2 + H 

5^4 X 10 11 

1/2 

0 

16. 

CHO + O = CO + OH 

5.4 x 10 U 

1/2 

0 

17. 

CHO + CH = CH + CO 
3 4 

2.5 x 10 10 

1/2 

0 

GO 

H 

CHO + OH = CO + HjO 

3 x 10 13 

0 

0 

19.* 

HCO + M = H + 00+M 

2 x 10 12 

1/2 

-14.4 x 10 3 


PLUS ALL REACTIONS IN TABLE II 


* These reactions are retained for purposes of 
comparing with t he mechanism of TABLE III 



TABLE V. - OPERATING CONDITIONS FOR TWO- PHASE 
RECIRCULATION ZONE COMBUSTION CALCULATIONS 


Fuel : 


C 8 H 18 


Injection Temperature Fuel: 300° K 

Air : 300°K 


Equivalence Ratio 
Fraction Liquid Fuel: 


1.0 

0.5 (for all cases 
except d 1 = 0) 


Pressure: 1 atm 

Variable: Droplet diameter and 

m/v (residence time) 
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TABLE VI 


VALUES OF 
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TABLE VII. - SWIRLING JET EXHAUST CHARACTERISTICS 


s 


0.112 

0.232 

0.06 

0.6 

u (msec) 


76.0 

103.0 

41.0 

52 .4 

o 

w (m/sec) 
mo 


16.7 

32.7 

4.8 

35.6 

fuel 


LPG 

LPG 

N/A 

N/A 

$ 


5.05 

5.05 

N/A 

N/A 

(a) . . 

fuel 


0.245 

0.245 

N/A 

N/A 

T ° K 


300 

300 

300 

300 

o 

►o' 

8 

1 

o* 

(kg / m 2 ) 

70 

19 

1 

81 

r . 

mm 

25 

25 

50 

50 

jet 

P m 

kg/m 3 

1 . 03 (10) 4 

1 . 03(10) 4 

1.03 (lCD 4 

1.03 (10) 4 

NOTE: For 
gas 

the combustion experiments, liquified 
(LPG) was used. 

petroleum 

The 

theoretical calculations used 

propane 

as the fuel 
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TABLE VIII. - ANNULAR JET INITIAL CONDITIONS 



Initial Data 


r//r wall 

U 

T 

°N2 

a 

02 

“fuel 


m/sec 

°K 


> 


.82 67 

61 

1500 

.768 

.232 

0.0 

.8309 

61 

1500 

.768 

.232 

o 

• 

o 

.8333 

61 

300 

0.0 

o 

• 

o 

1.0 

.8338 

61 

300 

0.0 

0.0 

1.0 

.8344 

61 

300 

o 

• 

o 

o 

• 

o 

1.0 

.8349 

61 

300 

o 

• 

o 

o 

• 

o 

1.0 

.8373 

61 

1500 

.768 

.232 

o 

• 

o 

.8414 

61 

1500 

.768 

.232 

0.0 


geometric constraints 

Case 1 p = cons tant = 10 atm.; outer wall contour variable, 
wail center body constant 


is variable 


outer wall 


= constant = 0 . 328m 
center body constant 
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FIGURE 1 GAS TURBINE COMBUSTOR FLOW SCHEMATIC 




FIGURE 2 - UNIFIED MODEL-MODULAR CONSTRUCTION 
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• - Calculated 


p = 1 atm 
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T 1 = 450° K 
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T 3 msec 
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FIGURE 5 - COMPARISON OF GASL THEORY WITH NOx EMISSIONS 
FROM A CO/AIR FIRED JET STIRRED REACTOR 
WITH H 2 ADDED IN THE AMOUNT OF ^ to 1% (VOLUME) OF CO 



p 


1 a tm 



With N + OH — — 
Reaction 



Without N+OH 
Reaction 



@ 


WSR @ T r 

Equi. @ T r 

Equi. /Adiab . 

WSR Incl.N+OH^NO+H 

Q WSR-Adiabatic (calculated) 
# WSR @ T r with 0.1% 

N in fuel (calculated) 

0 Experimental Data 
Ref. 12 


O 



Residence Tin 


FIGURE 7 - EARLY 



OF 0 







- Experimental (Ref. 25) 






RECIRCULATION ZONE ANALYSIS 
WITH FINITE-RATE CHEMISTRY 



Figure 10 - Effect of residence .time upon the fuel droplet 

diameter in the recirculation zone as a function 
of initial diameter 



Temperature 


RECIRCULATION ZONE ANALYSIS 
WITH FINITE- RATE CHEMISTRY 



Figure H - The effect of residence time upon the recirculation 
zone temperature as function of initial fuel droplet 
diameter . 
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PPM NO 


RECIRCULATION ZONE ANALYSIS 
WITH FINITE-RATE CHEMISTRY 



Initial Fuel Droplet Diameter, d 1 , (microns) 

Figure 12 - The effect of residence time upon the recirculation 
zone concentration of NO as a function of initial 
fuel droplet diameter 




% Water,* to Fuel (lb/lb) 


82 FIGURE 13- COMPARISON, PRESENT THEORY/TURBINE COMBUSTOR DATA - 

EFFECT OF WATER AND STEAM INJECTION UPON NO x REDUCTION 
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FIGURE 17 — COMPARISON. THEORY/EXPERIMENT, TANGENTIAL VELOCITY 
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FIGURE 19 


COMPARISON, THEORY/EXPERIMENT 








Swirl 

Factor = 0.232 

* Prediction, Finite Rate 
Chemistry 

Exp. Data, Chervinsky, 

Ref. 33 

See Table VII 
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5 . .... 6 7 


TEMPERATURE PROFILE 





FIGURE 23 - COMPARISON, THEORY/EXPERIMENT, 



TEMPERATURE VARIATION 
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FIGURE 24 - COMPARISON. THEORY/EXPERIMENT, AXIAL 


temperature 


VARIATION 
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FIGURE 26- TEMPERATURE PROFILE DEVELOPMENT WITHIN 
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FIGURE 27- TEMPERATURE PROFILE DEVELOPMENT WITHIN ANNULUS, CASE 2, r = CONSTANT = 0.328m 

wall 
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